Isotropic-nematic transition in liquid crystals confined between rough walls 



in 
o 
o 

(N 
> 

o 
oo 



o 

•4-* 

c3 



I 

c 

o 
o 



> 

00 
On 



in 
o 



c 

o 
o 



13 



D. L. Cheung 1 and F. Schmid 2 

1 Department of Physics and Centre for Scientific Computing, 
University of Warwick, Coventry, CV4 7AL, f/fQ 
2 Fakultat fur Physik, Universitat, Bielefeld, 33615 Bielefeld, Germany 

The effect of rough walls on the phase behaviour of a confined liquid crystal (LC) fluid is studied 
using constant pressure Monte Carlo simulations. The LC is modelled as a fluid of soft ellipsoidal 
molecules and the rough walls are represented as a hard wall with a number of molecules randomly 
embedded in them. It is found that the isotropic-nematic (IN) transition is shifted to higher pressures 
for rougher walls. 



I. INTRODUCTION 

The phase behaviour of many fluids alters when they 
are confined 1, . This is particularly true for ordered fluids 
such as liquid crystals (LC} 2 -. As well as being interesting 
in its own right, the effect of confinement LCs is of tech- 
nological significance. Many applications of LCs depend 
crucially on this interaction. 

Because of this interest there have been many stud- 
ies of confined LCs. Simulations have studied many as- 
pects of the effect of confinement on LCs, such as wetting 
behaviour— and anchoring strength^. The effect of con- 
finement on the phase behaviour has also been studied 5 - 6 . 
In a recent paper the effect of rough walls on the be- 
haviour of confined LCs was examined^. Specifically the 
structure of the fluid and the surface anchoring was stud- 
ied. It was found that the surface anchoring became 
weaker as the surfaces became rougher, which may have 
implications for the phase behaviour of the confined fluid. 
The aim of the present paper is to study the effect of 
rough surfaces on the isotropic-nematic (IN) transition 
of a confined LC. 



II. MODEL AND SIMULATION DETAILS 

In order to simulate both a large number of fluid 
molecules for a reasonable range of pressures and 
wall roughnesses simple models for both the molecule- 
molecule and molecule-wall interactions are used. 

A simplified version of the Gay-Berne (GB) potential 
is employed. The interaction between two molecules i 
and j, with positions and Vj, and orientations Uj and 
Uj is given by 

V(t ij ,u i ,u j ) = i u ^ q 



otherwise 



where ry = r s ; — rj and 
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per. cr(fy Ui, Uj) is the shape function given by^ 



a(rij,Ui,Uj) = a \l 



1 + XUi.Uj 



1 - xu 4 .Uj 
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r ij = \ r ij\i an d v%j = Tij/rij. (To is the molecular width, 
which defines the length scale used throughout this pa- 



where x = ~ 1)/( k2 + 1); with the elongation k = 3. 
This approximates the contact distance between two el- 
lipsoids, eo and (To define the energy and length scales 
used throughout this paper, with other quantities given 
in terms of these. Specifically the reduced pressure is 
P* = Pctq/cq and the reduced temperature is T* — 
UbT '/erj. As the potential is fully repulsive, temperature 
is not a significant variable and is set T* = 0.5 through- 
out. 

The walls are represented by a hard core potential act- 
ing on the centres of mass of the fluid molecules, giving 
rise to homeotropic anchoring. Roughness is introduced 
by embedding some molecules in the wallt These are 
given random positions in the xy plane and random ori- 
entations which are held fixed during the simulations. 
It should be noted that these are used as a convenient 
method to introduce disorder into the surface-molecule 
interaction and do not correspond to real molecules. The 
roughness of the wall is characterised by the density of 
these embedded molecules S. 

As the aim of this study is to investigate the phase 
behaviour, constant NPT simulations are more appro- 
priate than constant NVT simulations^ -. To retain a 
constant S, volume moves involve only changes in the z 
box length. The phase behaviour is probed using a se- 
ries of simulations at different pressures, starting from a 
high pressure P* = 2.50 in the nematic phase. The sim- 
ulated systems comprised 1200 fluid molecules and up to 
60 molecules in the walls. The fixed area of each wall was 
154.9 cr 2 . The cell width was in the range 22 - 28 <7q. To 
ensure some sampling of surface configurations three dif- 
ferent walls were studied for each value of S. The order 
parameters and directors for these are presented in Tab. 
□ 

Table 1. 

Order parameter and directors for the wall 
configurations studied. 
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Due to the number of different systems considered (3 
for each non-zero E and one smooth wall) only a coarse 
sampling of pressure is attempted; the pressure being 
changed in steps of 0.02 (in reduced units) near the NI 
transition. This limits the accuracy at which transition 
pressures may be determined. While techniques such 
as Gibbs Ensemble Monte Carlo (GEMC)±±*±£ or grand 
canonical Monte Carlo (GCMC) 10 simulations may be 
used to investigate phase equilibria the densities involved 
lead to unacceptably low acceptance rates for particle 
transfers. Studying more elongated moleculesi&ii (with 
lower pni) or softer interaction models^ would allow the 
use of these techniques. 

As the nematic phase has oricntational order which 
is lacking in the isotropic phase, the phase of the system 
may be determined by calculating the orientational order 
parameter. This is given by the largest eigenvalue of the 
ordering tensor 



Qa/3 = ^2 ( o U i" U «/3 ~ 9<W ) ( 4 ) 



where 5 a [3 is the Kronecker delta function. The order 
parameter S2 is then determined by diagonalising this. 
For a nematic S2 is typically > 0.4, while lower values 
are typical of an isotropic state. 

As the presence of surfaces gives rise to an inhomo- 
geneous distribution of molecules it is useful to consider 
quantities such as the density or order parameter in dif- 
ferent regions of the simulation cell. As previous studies 
indicate that a layer of well ordered molecules forms even 
for rough walls£, the order parameter in the central half 
of the cell S% ulk (l z /4 to 3/ 2 /4) may be taken to be the 
most useful variable for characterising the phase of the 
confined fluid The variation of density and order may be 
examined more closely by the density and order param- 
eter profiles. 




p ~ ~ A ' ' " p 



FIG. 1: Variation of bulk order parameter with pressure. 
Solid line with filled circles shows E = 0.0 (expansion), dotted 
line with open circles shows E = 0.0, dashed line with filled 
diamonds shows (a) E = 0.1, (b) E = 0.2, (c) E = 0.3, and 
(d) E = 0.4 (expansion) and the long dashed line with open 
diamonds shows (a) E = 0.1, (b) E = 0.2, (c) E = 0.3, and 
(d) E = 0.4 (compression). 



III. RESULTS AND DISCUSSION 
A. Phase behaviour 

Shown in Fig. ^ S% as a function of P. As can 
be seen, on expansion gradually decreases until 

P* = 2.25 (S\ ulk w 0.4), before falling to an isotropic 
state at P* = 2.20. For this system the bulk isotropic- 
nematic transition occurs at P N ibuik = 2.3Qi&. Note that 
due to finite size effects the order parameter does not go 
exactly to Q-'' '. The system shows a similar transition on 
compression for an isotropic starting state. The lower- 
ing of Pni relative to the bulk phase is caused by the 
formation of a highly ordered layer near the walls 6 ' 18 . 

This behaviour is similar to the rough wall with £ = 
0.1 (Fig. nji) . However, the NI transition seems to be 
shifted to a higher pressure (between P* = 2.25 and P* = 
2.30). Away from the transition the behaviour of is 
similar to that of the smooth wall system. For pressures 
above the PbuIuni the fluid in the cell bulk would be 
nematic anyway so the rough walls would have little effect 
on S bulk in this case, while for P much less than Pni 
the anchoring effect of the walls is insufficient to lead to 
ordering in the cell bulk. 

The effect of further increasing the surface roughness 
may be seen in Fig. which shows against P* 

for the £ = 0.2 wall. As may be seen the NI transition 
is shifted to a high pressure compared to the £ = 0.1 
wall (P NI « 2.30). Figure ^c) and (d) show the phase 
diagram for walls with E = 0.3 and E = 0.4 respectively. 
While Pni increases on going from E = 0.2 to E = 0.3, 
there is a much smaller difference between E = 0.3 and 
E = 0.4, suggesting that the effect of the rough walls 
becomes saturated with increasing roughness, similar to 



FIG. 2: (a) Density profiles for the smooth wall system. 
Solid line is P* = 2.25 (expansion), dotted line is P* = 2.20 
(expansion), dashed line is P* = 2.22 (compression), and dot- 
dashed line is P* — 2.26 (compression) (b) Order parameter 
profiles for smooth wall system. Symbols as in (a). 



FIG. 3: (a) Density profiles for rough wall with E = 0.2. 
Solid line is P* — 2.30 (expansion), dotted line is P* — 2.28 
(expansion), dashed line is P* = 2.28 (compression), and dot- 
dashed line is P* = 2.32 (compression) (b) Order parameter 
profiles for rough wall with E = 0.2. Symbols as in (a). 



the anchoring coefficient^. The change in order param- 
eter through the NI transition also appears to become 
smoother as £ increases. This may be due to averag- 
ing over different wall configurations as these may have 
different Pni- This change in the location of the IN tran- 
sition may also be seen from the large error bars for S\ ulk 
for pressures near Pni- 

The approximate transition pressures and order pa- 
rameters are presented in Tab. [21 These generally in- 
crease with surface roughness. This may be attributed 
to the rough walls disturbing the highly ordered surface 
layer, which disrupts the order in the bulk. This can 
also be seen from the decrease in the surface anchoring 
strength as £ increases^. 

Tabic 2. 

Approximate transition pressures and order parameters. 
Errors in the final decimal place are in parenthesises 



r) ex P Tj ex P Tjcmp ncmp nNexp 



Qlexp 
a 2 



^2 ^2 



0.0 2.25 2.20 2.26 2.22 0.41(4) 0.18(6) 0.40(7) 0.20(6) 

0.1 2.30 2.25 2.30 2.25 0.50(6) 0.15(7) 0.45(1) 0.15(7) 

0.2 2.35 2.25 2.34 2.30 0.57(1) 0.15(7) 0.52(9) 0.20(8) 

0.3 2.40 2.28 2.40 2.30 0.61(3) 0.14(5) 0.63(2) 0.16(2) 

0.4 2.35 2.28 2.40 2.28 0.52(7) 0.16(7) 0.62(9) 0.14(5) 

B. Fluid Structure 

Shown in Fig. [3^ are density profiles for the smooth 
wall system. These are shown for pressures just above 
and below the isotropic-nematic transition for both ex- 
pansion and compression (for the expansion runs these 
are P nem = 2.25 and Pi SO = 2.20 and for the compres- 
sion runs Pnem = 2.26 and Pi SO = 2.22). As can be seen 
there is a noticeable decrease in the height of the density 
peak at around z pa 2.5 when going from the nematic to 
isotropic phases (P* — 2.25 — > 2.20). However, on going 
from the isotropic to nematic phase (P* = 2.22 — > 2.26) 
there is a much smaller change in the density profile. 

The order parameter profiles for these pressures are 
shown in Fig. 03- For all pressures there is a well or- 
dered layer near the wall, which has similar values for 



all pressures. More variation is seen around the second 
peak z pa 2.5. For the isotropic fluids the order parame- 
ter profile is typically less than 0.1, while for the nematic 
S\ ulk (z) fa 0.4. 

Density and order parameter profiles for the rough wall 
with £ = 0.2 are shown in Fig. [3] The expansion runs are 
given at P nem — 2.30 and Pi SO — 2.28 and for the com- 
pression runs Pnem — 2.28 and Pi so = 2.32. Compared to 
the smooth wall the density profiles are less structured, 
with only a single peak at z ~ 2.5 visible for the low- 
est pressures shown, which also appears wider than the 
equivalent peak for the smooth wall, due to the disrup- 
tion of the layer structure due to the surface disorder—. 

The effect of the rough walls is more prominent for 
the order parameter profiles (fig. Near the wall 

(z)is significantly lower than for the smooth wall, 
due to the disordering effect of the wall molecules. The 
order minima is deeper (the minimum valuer 0.1 for all 
pressures shown) and is moved closer to the wall (z pa 
1.1). However, this does not extend into the bulk fluid. 

Density and order parameter profiles for £ = 0.4 are 
shown in Fig. 0] Profiles for the expansion runs are at 
Pnem = 2.35 and P* so — 2.28 and for the compression 
runs P* em = 2.28 and P* so = 2.40. The density profile 
shows even less structure than for £ = 0.2, with only 
a density peak at z pa 2.5 visible. However, this peak 
has a broad shoulder starting at z pa 2.0, corresponding 
to molecules lying in the plane of the wall in the region 
between the two layerst This is similar to the behaviour 
seen in simulations of smectic LCsSi. The effect of this 
bimodal distribution may be seen by the depletion of 
the order parameter in this region and by a peak in the 
biaxiality profile (c.f. Fig 2d of Refi). 



IV. CONCLUSIONS 

The phase behaviour of a model liquid crystal confined 
between two rough walls has been studied. Particular 
attention is paid to its behaviour in the vicinity of the 
isotropic-nematic transition. 

For smooth walls it is found that the nematic-isotropic 
transition is shifted to lower pressures compared to the 
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FIG. 4: (a) Density profiles for rough wall with E = 0.4. 
Solid line is P* = 2.35 (expansion), dotted line is P* = 2.28 
(expansion), dashed line is P* = 2.28 (compression), and dot- 
dashed line is P* — 2.40 (compression) (b) Order parameter 
profiles for rough wall with E = 0.2. Symbols as in (a). 

bulk. This may be attributed to the anchoring effect of 
the walls. As the walls become rougher, this effect is 
weakened so the transition is shifted to higher pressures. 
For pressures away from the transition pressure the rough 
walls appear to have a smaller effect. 

On examining the structure of the confined fluid, it 
is found that the density distribution changes little with 
pressure, rather changes in this are largely due to in- 
creasing wall roughness. The behaviour of the order pa- 



rameter near the wall is likewise determined by the wall 
roughness. 

To summarize, we find that surface roughness not only 
decreases the nematic order in the vicinity of the surface, 
but also shifts the ncmatic-isotropic transition to higher 
pressures in thin slabs. It would be interesting to study 
other quantities, in particular, the anchoring direction. 
On flat surfaces, the anchoring direction is homeotropic. 
On rough surfaces, we observe that the system sometimes 
switches to planar anchoring. However, this depends sen- 
sitively on the particular realization of the disorder, and 
much more extensive simulations of much larger systems 
would be necessary to clarify this issue. 
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